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We formulate a theory to the 2D Hubbard model in a framework free of finite size effect and 
numerical analytical continuation, yet containing the essential features of the 2D Hubbard model, 
i.e., the correct atomic limit for large frequencies and 2D spin fluctuations. As temperature is 
decreased for a 2D half-filled band, 2D critical fluctuations give rise to a strong local maximum in 
| Im'S(kF,ui) at u> — leading to a split peak in the spectral function. As U is increased, four peaks 
associated with antiferromagnetic and Hubbard bands begin to develop in small and intermediate 
frequency regimes. 

PACS numbers: 71.10.Fd, 71.27.+a 



Recently the Hubbard model has received consider- 
able attention, since it is believed to contain the essential 
physics of strong electron correlations found in the high 
temperature superconductors 0j . Although it was solved 
exactly in one-dimension J2| , the exact solution in higher 
dimensions is not known yet. For a half-filled 2D band, 
the T = ground state is believed to be antiferromag- 
netic insulating for all U, while at finite temperatures 
strong spin fluctuations destroy long-range antiferromag- 
netic order due to the Mermin- Wagner theorem ||. 

In the absence of a small expansion parameter of the 
model in the physically relevant regime, quantum Monte 
Carlo (QMC) simulations have played an important role 
in elucidating various dynamical properties at finite tem- 
peratures M. In spite of its exact nature, QMC calcu- 
lations suffer from two serious limitations, namely, small 
lattice size and numerical analytical continuation from 
imaginary frequencies to real ones. In fact, depend- 
ing on lattice size and on the method used to extract 
A(k,w) from G(k,iiv n ) computed by QMC, conflicting 
results have been reported for the single particle spectral 
function along the antiferromagnetic zone boundary for 
a half-filled 2D band with U = At and T = 0.1 - 0.25i 
HU ■ Although calculations based on the lowest order di- 
agram (second order in U) show the correct atomic limit 
for large iv, they fail to describe the correct 2D spin fluc- 
tuations. This atomic limit is reflected as the correct 
asymptotic behavior of the self-energy at large frequen- 
cies, J7 2 n/2(1 — n/2)/iu> n (in the paramagnetic state), 
which is necessary to produce the Hubbard bands at large 
enough U. Recently fluctuation exchange (FLEX) ap- 
proximation was applied to the 2D Hubbard model by 
Bickers et al. J?J and later extensively used by Bickers 
et al. H and by Dahm et al. 0. Despite its successful 
application to the D-wave superconductivity and many 
other issues, it was pointed out by Vilk and Tremblay 
pOf that FLEX approximation does not show both the 



precursors of antiferromagnetic and Hubbard bands re- 
sulting from the correct 2D spin fluctuations and atomic 
limit for large u, respectively. This is because FLEX 
approximation completely neglects the frequency depen- 
dent vertex corrections for the self-energy within its self- 
consistent structure. While the 2D spin fluctuations were 
quite successfully described by Vilk and Tremblay by 
imposing the exact sumrules on the spin and charge sus- 
ceptibilities, their self-energy does not show the correct 
asymptotic behavior at large frequencies. Therefore, it is 
highly desirable to formulate an approximation scheme 
for the 2D Hubbard model in a manner free of finite size 
effect and numerical analytical continuation, yet contain- 
ing the essential features of the 2D Hubbard model, i.e., 
the correct atomic limit for large lo and 2D spin fluctu- 
ations. In this Letter, we propose such a theory for the 
first time. 

We begin with reconsideration of paramagnon theory 
first proposed by Berk and Schrieffer fyj]. This theory 
has some good features, namely, explicitly enforcing the 
rotational invariance of the spin susceptibilities together 
with including the particle-particle channel, and the cor- 
rect prediction of the zero temperature phase transition 
for any U in the half-filled 2D band. This theory, how- 
ever, always gives rise to a finite temperature antiferro- 
magnetic instability for large enough U due to the insuf- 
ficient treatment of strong 2D spin fluctuations, in con- 
flict with the Mermin- Wagner theorem. At frequencies 
larger than the bandwidth W = 8, the self-energy ob- 
tained from this theory is not guaranteed to give the 
correct asymptotic behavior. Recently we found that 
within the same structure for the self-energy and var- 
ious susceptibilities as that in the paramagnon theory, 
both the Mermin- Wagner theorem and correct atomic 
limit for large 10 can be satisfied simultaneously by intro- 
ducing renormalized interaction constants U sp , U c h, and 
U pp in the spin, charge, and particle-particle channels, 
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respectively. 

In our formulation the self-energy is expressed in terms 
of these renormalized interaction strengths to be deter- 
mined later: 
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T, N and U are the absolute temperature, number of 
lattice sites and Coulomb repulsion energy, k is a com- 
pact notation for (k, iuj n ) where ioj n are either Fermionic 
or Bosonic Matsubara frequencies. G°(k) is the noninter- 
acting Green's function and Xph(q)i X PP (q) are irreducible 
particle-hole and particle-particle susceptibilities, respec- 
tively, which are computed from 
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Dynamical spin, charge and particle-particle susceptibil- 
ities are calculated by 
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We determine U sp , U c h, and U pp by imposing the follow- 
ing three exact sumrules to Eq. (H): 
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where the Pauli exclusion principle {n 2 ) = (n a ) is ex- 
plicitly used in the spin and charge channels. By finding 
U sp , U c h and U pp through Eq. (||), any possible magnetic 
instability can happen only at zero temperature, because 
the right hand sides are always finite, in consistency with 
the Mermin- Wagner theorem. The first two sumrules in 
Eq. U were previously used by Vilk and Tremblay 



(b) 



0.20 - 



0.10 - 



(1) 



0.00 



-8.0 



0.0 
CO 



8.0 



0.00 




FIG. 1. (a) Spectral functions at k = (iv/2,n/2) and (b) 
density of states, for (7 = 4 at T =0.10, 0.15 and 0.20 denoted 
as the solid, dotted and dashed curves, respectively. 



to study collective spin and charge fluctuations in the 
2D Hubbard model. The correct asymptotic behavior of 
the self-energy at large frequencies can be easily checked 
by substituting Eq. (^) and {n^n^) — n 2 /4 for nonin- 
teracting electrons into Eq. p). By making an ansatz 



U sp = U(n^n^)/({nij)(ni)) |12j, the double occupancy is 
computed self-consistently in the first sumrulc in Eq. (Q) 
and U c h, Upp can be obtained by the other two sumrules. 
Throughout the calculations the unit of energy is t. We 
used a 128 x 128 lattice in momentum space and per- 
formed the calculations by means of well-established fast 
Fourier transforms (FFT). It should be also noted that 
we used a real frequency formulation in Eqs. (0)-(|J) to 
avoid any possible uncertainties associated with numeri- 
cal analytical continuation. 

As a first application of the present theory to the 
2D Hubbard model, we study the half-filled 2D Hub- 
bard model in the intermediate coupling regime (U = 4) 
where QMC calculations have shown conflicting results 
for the single particle spectral function at the noninter- 
acting Fermi surface [g|^]. In this Letter we would like 
to give a definitive answer to this issue for the first time. 
The spectral functions at k — (ir/2, ir/2) and the density 
of states are presented in Fig. |l|(a) and (b) for U = 4 
at T = 0.1, 0.15, and 0.2. At T = 0.2 (dashed curve 
in Fig. |l](a)) a single quasiparticle peak is found at the 
Fermi energy. As the temperature is slightly decreased 
to T = 0.15 (dotted curve), the single particle peak be- 
gins to split into two, leading to a pseudogap Jll| at the 
Fermi energy. This is a precursor of antiferromagnetic 
bands resulting from strong 2D spin fluctuations at low 
temperatures. At T = 0.1 (solid curve), the two peak 
structure with a pseudogap inside becomes much more 
pronounced. The present results support the existence 
of a pseudogap in the spectral function at T < 0.2 for 
U = 4. For all three temperatures, however, the density 
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FIG. 2. (a) Imaginary and (b) real parts of the self-energy 
at k = (7r/2,7r/2) for U = 4 at T =0.1, 0.15 and 0.2 denoted 
as the solid, dotted and dashed curves, respectively. The solid 
straight line in (b) is u> — e(k). 

of states (Fig. 0(b)) shows a pseudogap at the Fermi en- 
ergy. In spite of a single peak structure in the spectral 
function for T = 0.2 (dashed curve), the density of states 
is significantly suppressed at the Fermi energy compared 
with that for non-interacting electrons which shows a log- 
arithmic divergence. Because of its accumulative nature, 
the density of states appears more sensitive to the change 
of a quasiparticle state at the Fermi energy than the spec- 
tral function itself does, at least, for a half-filled 2D band. 
This situation for T = 0.2 is the best result found by 
Dahm et al. || within FLEX approximation in an effort 
to see the precursor effect of the antiferromagnetic bands. 
The small bumps at ±(4 — 5)i both in the spectral func- 
tion and density of states are precursors of the Hubbard 
bands, which will be discussed later. 

This anomalous behavior is understood more clearly by 
examining the imaginary and real parts of the self-energy 
shown in Fig. |](a) and Fig. ||(b) , respectively for the same 
parameters as in Fig. [ij At T = 0.2 (dashed curve), 
the imaginary part of the self-energy becomes smaller in 
magnitude as the Fermi energy is approached, indicat- 
ing a Fermi liquid-like behavior. As the temperature is 
further decreased to T = 0.15 (dotted curve) to T = 0.1 
(solid curve), however, a drastic change happens near 
the Fermi energy. Due to 2D critical fluctuations, the 
scattering rates at the Fermi energy grow exponentially 
as ~ £ ~ exp(constant/T). This singular scattering by 
exchange of strong 2D spin fluctuations is responsible 
for the strong suppression of the spectral weight at the 
Fermi energy, leading to the formation of a pseudogap in 
the spectral function. 

The corresponding real part of the self-energy obtained 
from the imaginary part by means of the Kramers-Kronig 
relations is shown in Fig. |^(b) . The peak condition in the 
spectral function is determined by the intersection of the 



FIG. 3. (a) Spectral functions at k = (n/2,ir/2) and (b) 
density of states, for U =4, 6 and 8 at T = 0.1 denoted as the 
solid, dotted and dashed curves, respectively. 



real part of the self-energy and w = e(k) denoted as the 
solid straight line in Fig. ^(b). e(k) = —2t(cosk x +cosk y ) 
for nearest neighbor hopping. At T = 0.2 (dashed curve), 
the peak condition is satisfied only at the Fermi energy 
with its slope being negative, a characteristic feature of 
the Fermi liquid. As the temperature is decreased to 
T = 0.15 (dotted curve) to T = 0.1 (solid curve), its 
slope becomes positive and larger than unity at the Fermi 
energy, and thus the peak condition is satisfied at three 
different locations. Because of large scattering rates at 
the Fermi energy, however, only two peaks appear in the 
spectral function, consistent with the results in Fig. ^. 

Since the present formulation for the 2D Hubbard 
model can describe quite reasonably both 2D critical spin 
fluctuations and the correct atomic limit for large ui, it is 
of great interest to examine how the precursors of antifer- 
romagnetic and Hubbard bands evolve as the interaction 
strength is increased. For U = 4 (solid curve in Fig. ||(a)), 
a split peak with developing small bumps at ±(4 — 5)i 
is found. As U is increased to 6 (dashed curve), the 
antiferromagnetic bands become significantly suppressed 
and the Hubbard bands grow substantially near ±6. For 
U = 8 (solid curve), the spectral weight inside the Hub- 
bard bands becomes further suppressed due to the large 
Coulomb repulsion, and as a result the spectral weight 
associated with the Hubbard bands becomes dominating 
over that with the antiferromagnetic bands. The Hub- 
bard bands for U = 8 occur at larger frequencies than 
±U/2, since the asymptotic behavior of the self-energy, 
U 2 n/2(\ — n/2)/ito n , sets in at much higher frequencies 
than the bandwidth. This is because the strong max- 
imum of the scattering rates at the Fermi energy sig- 
nificantly enhances the real part up to high frequencies. 
The appearance of four peaks in the spectral function for 
U = 8 is consistent with recent QMC calculations 15 
except a difference in the relative strength of the two dif- 
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FIG. 4. (a) Imaginary and (b) real parts of the self-energy 
at k = (7r/2,7r/2) for U =4, 6 and 8 at T = 0.1 denoted as 
the solid, dotted and dashed curves, respectively. The solid 
straight line in (b) is iv — e(k). 



ferent bands. This feature is also similar to the one from 
the 1/d effect in the infinite dimensional Falicov-Kimball 
model |L6| . The density of states obtained from the av- 
erage of the spectral function in the whole Brillouin zone 
shows more suppressed antiferromagnetic bands and thus 
more pronounced Hubbard bands in Fig. ||(b) than the 
spectral function does at (ir/2,ir/2). 

As the interaction strength U is increased, the imag- 
inary part of the self-energy grows rapidly both at the 
Fermi energy and in the intermediate frequency region 
(~ ±4) shown in Fig. ^(a). The former comes from 
strong 2D critical fluctuations near half-filling and the 
latter mainly from the lowest order diagram (second or- 
der in U) which is responsible for the Hubbard bands for 
large enough U. As U is increased to 8 (dashed curve in 
Fig. (l|(b)), seven solutions are found in the peak condi- 
tion. The most outer two solutions lead to the Hubbard 
bands and the most inner two near the Fermi energy 
to the antiferromagnetic bands, and the other three ac- 
companied by large scattering rates yield the incoherent 
background. This result is similar to the earlier report by 
Kampf and Schrieffer [H| of the developing multiple so- 
lutions in the peak condition by increasing the spin-spin 
correlation length. Besides the total number of solutions, 
however, there are some qualitative differences between 
these two results. Because of the insufficient treatment 
of the 2D critical fluctuations in their phenomenological 
antiferromagnetic spin fluctuation spectrum, their solu- 
tion at the Fermi energy shows a single quasiparticle 
peak, while in our calculations a quasiparticle state is 
destroyed for a half-filled band. Due to the same reason, 
the pseudogap found in their paper does not come from a 
suppressed spectral weight inside the antiferromagnetic 
bands but instead inside the Hubbard bands. 

In summary, we formulated a theory to the 2D Hub- 



bard model in a manner free of finite size effect and nu- 
merical analytical continuation, yet containing the essen- 
tial features of the 2D Hubbard model, i.e., the correct 
atomic limit for large u> and 2D spin fluctuations. As 
the temperature is decreased for a 2D half-filled band, 
anomalous behaviors are found near the Fermi energy 
such as a split peak in the spectral function, a large 
positive slope greater than unity in the real part of the 
self-energy and a strong local maximum in the scattering 
rates. As the interaction strength U is increased, four 
peaks associated with the antiferromagnetic and Hub- 
bard bands begin to develop in the small and intermedi- 
ate frequency regimes. 
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